Shift of BEC Temperature of Homogenous Weakly Interacting Bose Gas* 
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We report on the computation of the shift of the Bose-Einstein condensation temperature for a 
homogenous weakly interacting Bose gas in leading order in the diluteness parameter an 1 / 3 , where 
a is the scattering length and n is the particle density. The perturbative series, which is afflicted by 
infrared divergences, is resummed by means of variational perturbation theory. Using coefficients 
through seven loops, we arrive at AT C /T C = 1.27±0.11 an 1 / 3 , which compares favorably with recent 
Monte-Carlo data. 
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I. INTRODUCTION 

The Bose-Einstein condensation (BEC) temperature 
T c of an ideal gas of spin-0 Bosons is given by 



T 



2tt 



C(3/2) 



2/3 



(1) 



where m is the mass of the bosons, n their number den- 
sity and we work throughout in units where fee = ft = 1. 
Infrared (IR) divergences prevent the perturbative evalu- 
ation of the shift of the BEC temperature due to a small 
repulsive interaction, parameterized by the correspond- 
ing s-wave scattering length a. Physically, these IR di- 
vergences correspond to critical fluctuations that change 
the universality class from Gaussian to that of a three- 
dimensional scalar 0(2)-symmetric field theory. 

Although it is clear that T c changes by a small amount 
if the dimensionless diluteness parameter an 1 / 3 is small 
(for small momenta, a is the only relevant parameter of 
the two-body interaction potential) , it was for a long time 
unclear, what the leading power in terms of this param- 
eter or the coefficient in front of it would be. Many at- 
tempts on the problem have provided different powers 
and coefficients [1-24]. Recently, however, it has been 
shown that the leading and next-to-leading behavior is 
given by [7, 8, 15] 



= dan 1 / 3 + [c' 2 lnfan 1 / 3 ) + c£] (an 1 / 3 ) 2 + • • • (2) 



In the work presented here, we determine c\ numerically 
(while c' 2 is known exactly, c? 2 ' may be determined by a 
non-perturbative calculation analogous to the one pre- 
sented here for determining ci; for the exact value of c' 2 



and a MC estimate of c 2 ', see Ref. [16]). Results for c\ 
in the literature range from —0.93 [9] to 4.7 [2], see Refs. 
[14, 17] and the review [25], also Fig. 1 below. In [7] it was 
shown that c\ is exclusively generated by infrared fluc- 
tuations and that consequently three-dimensional field 
theory is sufficient for the determination of ci. It ap- 
pears that the most reliable results for c± before our work 
have been obtained by Monte-Carlo (MC) simulations in 
three-dimensional field theory [13, 14]. 

Here we expand AT C in a perturbative series, which 
is subsequently resummed to obtain the physical non- 
perturbative limit. While simple schemes such as Pade 
approximants may slowly converge towards the correct 
value, our scheme of choice is Kleinert's field theoretic 
variational perturbation theory (VPT, see [26-28] and 
Chapters 5 and 19 of the textbooks [29] and [30], respec- 
tively; improving perturbation theory (PT) by a vari- 
ational principle goes back at least to [31]), which has 
proven very useful for the investigation of critical phe- 
nomena. 



II. FIELD THEORY 

The thermodynamic equilibrium properties of a ho- 
mogenous gas of spin-0 Bosons may be described in 
the grand canonical ensemble with the help of a non- 
relativistic (3 + l)-dimensional field theory in imaginary 
time r, given by the Euclidean action 
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(3) 

where /1 is the chemical potential. This assumes that the 
momenta of the particles are small enough so that their 
interaction potential is well described by just one pa- 
rameter, the s-wave scattering length a. It also assumes 
that three- and more particle interactions are rare, i.e. 
that the gas is dilute. It turns out that the leading per- 
turbative contribution to AT C is oc a 2 . In consequence, 
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the leading contribution to AT C arises exclusively from 
terms that are infrared divergent in the framework of 
PT. This is why, for the determination of c\, we can 
set all Matsubara frequencies to zero or, equivalently, 
work with a three-dimensional field theory [7]. Denot- 
ing the zcro-Matsubara modes by ipo, we define the fields 
and parameters of this theory by ipo — ^/ m T{4>\ + %4> < i), 
^baro = —2mfi, u — AS-KamT, and obtain the three- 
dimensional Euclidean action 
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For the determination of quantities that are not governed 
by the leading IR divergences, one may still use a three- 
dimensional field theory which, however, is obtained by a 
more complicated matching procedure, which introduces 
corrections to the above relations of the parameters of 
the (3 + l)-dimensional and the three-dimensional theory 
and also necessitates the inclusion of more interaction 
terms. This matching is described in detail in [16]. 

From the ideal-gas result (1) one obtains, in leading 
order, the relation [7] 



T 



2 An 
'3~rT' 



(5) 



where AT C is the shift of the condensation temperature 
for fixed n and An is the shift of the critical particle 
density for fixed condensation temperature. In our field- 
theoretic setup, An is given by 



An = A(ip*ip) = mTA(<j) 2 }. 
Combining (1), (2), (5) and (6), ci is given by 

A(0 2 ) 



ci 



a 



Nu 



(6) 



(7) 



with a = -256tt 3 /[C(3/2)] 4/3 w -2206.19 and the re- 
maining task is to compute the critical limit of 



2 > 



Nu 
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[G(p)-G (p)] 
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(8) 

where G and Go are the interacting and the free propa- 
gator, respectively. In (7) and (8) we have conveniently 
generalized the model (4) to an O(N) field theory, where 
4> = (0i, . . . , 4>n)j 4 }2 = 0a0a- This allows to make con- 
tact with the exactly known large- N result for c\ [8] and 
with the MC results for N = 1 and N = 4 [32]. 

The theory (4) is superrenormalizable. A convenient 
and popular renormalization scheme is given by defin- 
ing a renormalized parameter r — rb ar c — ^(0), where 
£(p) is the self-energy. Now all Feynman diagrams can 
be computed without rcgularization. The critical limit 
is obtained by letting r — > 0, as can be seen from 
the structure of the full propagator, which reads now 
G(p) = l/{p 2 +r- [E(p) - S(0)]}. For the free propa- 
gator follows G Q (p) = l/{p 2 + r) [20]. 



In a scheme used in [12, 18, 19, 22, 33], the free prop- 
agator is taken always critical, Go(p) = 1/p 2 - This gen- 
erates an unnatural non-zero one-loop contribution to 
A(0 2 ), even in the absence of interactions. Although 
this contribution vanishes as r — > 0, it strongly influ- 
ences resummation and is responsible for the small value 
of ci = 0.92 ± 0.13 [34] in Ref. [22]. 



III. PERTURB ATIVE SERIES FOR a 



In PT, A{<j) 2 )/Nu may be written as a power series 
in u r = Nu/Airr 1 / 2 . 
ci = lim^^oo ci(u r ), 



Defining a function c\{u r ) with 
we have 



ci{u r ) = a'^2 a l u 'r 



(9) 



i=i 



The inclusion of N in the definition of u r is motivated by 
the fact that then the ai remain finite even in the limit 
N — > oo, facilitating the comparison with the exactly 
known large- TV result [8]. The successful application of 
VPT to the large- N limit has been demonstrated in [23] . 

The perturbative series can be represented by Feyn- 
man diagrams and then becomes a loop expansion. The 
subtraction of the zero- momentum part of the self-energy 
in the full propagator has to be performed recursively for 
all subdiagrams. We denote this procedure by an opera- 
tor TZ [23, 24] . The first non-zero perturbative coefficient 
arises at the three-loop level, since the lowest momentum- 
dependent self-energy contribution has two loops. The 
perturbative coefficient at L loops is given by 



a L u 



L-2 
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(10) 



where til is the number of contributing diagrams in L- 
loop order and D^-k is the fc-th L-loop diagram. There 
are Ul — 0,0,1,1,5,12,56 diagrams at one- through 
seven-loop order (see [35] for a convenient way of con- 
structing all necessary diagrams together with their 
weights) and the expansion for c\(u r ) starts out as 



ci(u r ) 



a 
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E.g., the coefficient 03 is given by 
a^Ur 



(11) 
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The results for diagrams through six loops in the renor- 
malization scheme employed here can be found in [36]. 
The seven-loop diagrams have been computed in 1991 
and used for the determination of critical exponents for 
TV = 0, 1, 2, 3 in [37]. They were provided to the author 
by B. Nickel [38] and may be found in [24]. The resulting 
perturbative coefficients for TV = 2 are a 3 = —1.01209 x 
10~ 4 , a 4 = 2.99626 x 10~ 5 , a, = -1.19872 x 10~ 5 , 



a 6 = 5.85519 x 10~ 6 , a 7 = -3.30467 x 10" 



IV. RESUMMATION AND RESULTS 

For obtaining the u r — > oo limit of c\(u r ), the series 
(9) needs to be resummed. For critical fluctuations phe- 
nomena, VPT is an appropriate tool, since it is designed 
to accommodate the large-w r behavior of physical quan- 
tities in the critical limit, given in our case by 



ci(u r ) = a y~] f„ 



(13) 



with an irrational exponent u>' . The leading class of cor- 
rections of a physical quantity that remains finite in the 
critical limit are integer powers of t MV [30, 39, 40], where 
t = (T — T c )/T c , v is the critical exponent of the corre- 
lation length and u = f3'(g*) in a renormalization group 
approach. In our renormalization scheme, the propagator 
obeys G(p = 0) = 1/r oc i~ 7 . Using the universal scaling 
relation 7 = v(2 — 77), where r\ is the anomalous dimen- 
sion of the critical propagator, i.e. G(r = 0) oc l/p 2 ~ v in 
the small-p limit, we see that we have to set 



2uj 
2~^ 



(14) 



in (13). 

The ansatz (13) docs not account for so-called con- 
fluent singularities which cause the true large-w r expan- 
sion to also contain other negative powers of u r , which 
are subleading at least compared to . We can ex- 
pect methods that can accommodate the leading behav- 
ior in (13) correctly to converge faster to the true result 
than methods having the wrong leading behavior such as 
e.g. Pade approximants or the linear S expansion (LDE), 
the latter being used extensively for the current prob- 
lem [12, 18-21]. On the other hand, convergence will be 
slowed by the fact that we do not make an ansatz reflect- 
ing the full power structure in u r , but the expansion (13) 
will try to mimic the neglected subleading powers. 

The alternating signs of the cil suggest that the pertur- 
bative series for c\(u r ) is Borel summable. In the context 



of critical phenomena, such series have been successfully 
resummed using Klcincrt's VPT. Accurate critical expo- 
nents [27, 28, 30] and amplitude ratios [41] have been 
obtained. For a truncated partial sum a i ul r ~ 2 °f 

(9), the method requires replacing 



,1-2 



{tu) l ~ 



1 + t 



-(*-2)/u/ 



(15) 



(note that this is an identity for t = 1), reexpanding 
the resulting expression in t through t L ~ 2 , setting t = 1 
and then optimizing in u, where optimizing is done in 
accordance with the principle of minimal sensitivity [42] 
and in practice means finding appropriate stationary or 
turning points. That is, we replace 



k=0 



-{I -2)1 J 
k 



- 1 



(16) 



and optimize the resulting expression in u. For u r — > 00, 
we obtain the L-loop approximation of /o, 
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opt fl 



L-l 



fe=0 



-{l-2)/u 
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(17) 

Results are only available starting at four loops, since 
two non-zero perturbative coefficients are necessary for 
VPT to work. E.g., in four-loop order the optimization 
yields 



/o 4) ~ °Ptfi 



03 1 + 



U + CLiU 



0§ 



4a4 



(18) 

as the best attempt at determining /o and thus c\ = afo- 
For N = 2, we have uj = 0.79 ± 0.01, 77 = 0.037 ± 
0.003 (see, e.g., [30]) and thus, using (14), J = 0.805 ± 
0.011. The four- through seven-loop results arc ci = 
0.948, 1.062, 1.126, 1.161. 

u>' may also be determined self-consistently [27, 29, 30]. 
Assuming a behavior of c\{u r ) as in (13), the quan- 
tity d\nc\(u r ) / dhiu r has an expansion of the same 
type with the same u>' as c\(u r ), but with a vanishing 
large-u r limit (i.e. its large-u r expansion starts out with 
/o = 0). ui' is tuned such that the value VPT gives for 
d\nci(u r ) / d\nu r is zero in a given loop order. This u)' 
is then used as an input for the determination of the ap- 
proximation of ci(u r ) at the same loop order. The five- 
through seven-loop results are c\ = 1.399, 1.383, 1.376. 

From the fact that the fixed-a/ results monotoni- 
cally rise with the loop order, while the results with 
self-consistent ui' fall, we estimate that the true value 
lies inbetween and conclude from the resummed results 
through seven loops that 



ci = 1.27±0.11. 



(19) 



Our treatment for N = 2 may easily be repeated for 
arbitrary N. E.g., for N — 1 one obtains 



ci = 1.07 ±0.10, 



(20) 
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while for TV = 4, the result is 

ci = 1.54 ±0.11. (21) 

V. DISCUSSION 

Many different results have been obtained for c\; see 
[24] and corresponding references therein for criticisms of 
most of them. A comparison of our results through seven 
loops with most other results found in the literature is 
given in Fig. 1. Our result (19) is in agreememt with the 
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FIG. 1: Comparison of c\ from VPT as a function of the num- 
ber of loops L for N = 2 with most results from other sources. 
Upper dots: self-consistent (sc) to' . Lower dots: lo' — 0.805. 
LDE* indicates use of LDE with a resummation method for 
accelerated convergence. The label VPT* indicates the inclu- 
sion of a non-zero one-loop term in Ref. [22], which is absent 
in the present treatment, labeled by VPT. 

apparently most reliable other sources available, namely 
the MC results from the three-dimensional theory: c\ = 
1.29 ±0.05 by Kashurnikov, Prokof'cv and Svistunov [13] 
and ci = 1.32 ± 0.02 by Arnold and Moore [14]. 

MC results for N = 1 and N — 4 are also available. 
X. Sun has computed A(<T) 2 ) c rit./w for these cases [32] 



which translate to c\ — 1.09 ± 0.09 for TV = 1 and 
ci = 1.59 ± 0.10 for N = 4. Our values (20) and (21) 
agree well with these results. This provides additional 
confidence that in fact both MC and our VPT results 
are trustworthy estimates for c\. 

We would like to emphasize that the coefficient c\ is of 
physical relevance not only for the case of a strictly ho- 
mogenous gas, but also for a trap that is sufficiently wide 
to allow for critical fluctuations [10, 43]. If expressed as a 
function of the total number of particles N p , the leading 
effect on the shift of the condensation temperature is ob- 
tained from mean field theory and only in second order 
the influcence of critical fluctuations can be felt [43, 44]. 
In contrast, if expressed in terms of the central density 
n of particles, the leading effect is determined by critical 
fluctuations as we have seen here. It would therefore be 
desirable to know not only the total number of particles 
N p (as in a recent experiment [45]) but also the central 
density of particles in experiments that measure the con- 
densation temperature. In this context we would also 
like to stress that it is no contradiction that the sign of 
the change of T c is different in both cases, since different 
physical quantities N p and n are kept fixed, respectively 
Let us turn this around: For a given condensation tem- 
perature, a small repulsive interaction causes the corre- 
sponding central density in a wide trap to be lowered 
and the total number of particles to be raised compared 
to the ideal gas case. 
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